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SUMMARY 

A  seven-month  study  was  conducted  with  $9^,000  of  19&i  Augmented  Independent 
Development  Pund^  to  determine  the  fundamental  nature  of  the  gas  dynamic  re¬ 
circulation  process  which  gives  rise  to  staging  force,  moment,  and  environment 
problems  in  nearly  all  rocket  powered  separations.  Parametric  cold  gas  experi¬ 
ments  were  conducted  for  a  series  of  progressively  more  complex  idealized  inter¬ 
stage  configurations  which  would  exhibit,  in  sequence,  the  steps  in  the  recircu¬ 
lation  phenomenon.  Theoretical  work  was  initiated  which  is  expected  to  lead  to 
a  method  for  waking  preliminary  design  estimates  of  recirculation  effects  on  a 
wide  variety  of  potential  vehicle  staging  configurations# 

The  experimental  program  was  completed  in  full,  with  sufficient  force,  pressure, 
and  photographic  data  taken  on  all  the  planned  configurations  to  support  theoreti¬ 
cal  analysis.  Plotted  pressure  and  force  data  taken  on  the  majority  of  the  experi¬ 
mental  configurations  is  presented  in  an  Appendix.  Seme  of  the  data  is  not  included 
due  to  lack  of  time  for  reduction  and  plotting.  A  comparison  between  certain  of  the 
data  and  previous  experimental  results  is  made,  showing  good  agreement. 

t 

A  discussion  of  theoretical  concepts  concerning  interstage  gas  dynamics  is  given, 
in  which  a  scheme  for  calculating  recirculation  is  envisioned  as  an  iteration  pro¬ 
cedure,  involving  a  series  of  independent  stagnations  and  intervening  supersonic 
re-accelerations  within  the  interstage.  Examination  of  the  data  shows  such  a  step¬ 
wise  process  of  stagnations  to  exist,  and  thus  the  envisioned  procedure  appears 
feasible,  providing  each  of  the  independent  steps  can  be  correctly  calculated. 

Theoretical  analysis  was  couple  ted  on  the  first  stagnation  process  in  the  sequence 
of  steps  leading  to  recirculation,  and  a  preliminary  method  for  the  second  stagnation 
within  the  interstage  was  devised.  However,  the  first  stagnation  method  was  found 
to  be  inadequate  within  a  high  pressure  interstage  environment,  when  it  was  compared 
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to  the  experimental  data  toward  the  end  of  the  project.  The  remaining  time  was 
devoted  to  studying  how  this  method  could  be  corrected  and  what  implications  could 
be  drawn  concerning  the  interstage  calculation  as  a  whole#  The  inadequacies  of 
the  theory  in  the  interstage  environment  are  examined  by  comparing  various  modi¬ 
fications  of  the  theory  to  the  experimental  data.  Ths  method  of  predicting  the 
separation  rocket  jet  boundary  appears  inadequate,  with  no  better  method  available 
which  is  also  sufficiently  rapid  for  the  computerized  interstage  application.  The 
assumption  of  zero  mass  and  momentum  flux  in  the  jet  boundary  region  appears  poor# 

While  the  proper  value  can  easily  be  assigned  to  jet  boundary  strength,  the  inter¬ 
action  of  the  jet  boundary  of  this  strength,  with  flow  attempting  to  leave  the  first 
stagnation  region  appears  to  invalidate  the  boundary  conditions  for  the  stagnation 
assumed  in  the  theory#  However,  under  experimental  conditions  where  the  jet  boundary 
is  very  weak,  the  first  stagnation  theory  does  shew  reasonable  agreement  with  experi¬ 
ment,  implying  that  the  basic  equations  of  the  theory  are  valid# 

From  the  comparisons  of  theory  and  experiment  it  is  found  that  the  interaction  of  the 
jet  boundary  with  the  first  stagnation  flow  field  produces  a  thin  gas  layer  along  the 
floor  of  the  interstage#  Thus  it  is  concluded  that  interstage  vents  located  at  floor 
level  would  be  particularly  effective  in  reducing  recirculation  for  single-stage  rocket 
powered  Reparation  designs#  Likewise,  sealing  the  region  of  each  interstage  floor,  in 
a  design  involving  simultaneous  rocket  separation  of  several  adjacent  stages,  may  be 
expected  to  greatly  reduce  recirculation  between  the  stages# 
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Croaa-aectlonal  area  (area  in  x-r  plane)  -  radii 

Height  of  pressure  tap  on  experimental  splitter  fence 
from  interstage  floor  plate  -  radii 

Nozzle  mass  or  momentum  constant  (see  Appendix  I) 

Height  of  pressure  top  on  experimental  circular  fence 
from  bottom  of  fence  -  radii 

Radial  distance  from  experimental  nozzle  centerline  to 
splitter  fence,  measured  normal  to  fence  -  radii 

Height  of  pitot  probe  from  interstage  floor  plate  -  radii 

Height  of  bottom  of  experimental  circular,  fence  from  inter¬ 
stage  floor  plate  -  radii 

Stagnation  layer  area  error  sensitivity  constant  (see  Appendix 

I) 

Mach  number  of  gas  flow 

Sides  of  wedge-shaped  control  volume,  V* 

Flux  of  mass  -  slugs/ sec* 

Pressure  -  lbs/in^ 

Radial  distance  from  nozzle-plate  centerline  to  an  experimental 
pressure  tap  -  radii 

Radial  distance  from  nozzle  centerline  -  radii 


s  A  point  on  the  interstage  floor  shock 

T  Thrust  parameter  -  lb. 

V  First  stagnation  control  volume,  bounded  by  the  floor  shock,  the 

wedge  sides  (m).  the  interstage  floor,  and  the  layer  exit  (d) 

-  see  Figure  (8) 

Volume  occupied  by  jet  flow  between  nozzle  and  interstage  floor 
shock  -  see  Figure  (8) 

x  Axial  distance  from  nozzle  exit  -  radii 
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Flow  angle  measured  from  jet  axis  -  degrees 
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_ _ Explanation 

Angle  of  tilt  of  experimental  nozzle-jet  axis  from 
plate  -  container  axis  -  degrees 

Separation  rocket  exhaust  ratio  of  specific  heats 

Distance  of  standoff  from  the  interstage  floor  surface 
radii 

Width  of  first  stagnation  layer  measured  in  the  x  direction  - 
radii 

Prandtl  -  Meyer  angle  -  degrees 

Angle  between  pitot  probe  axis  and  the  plane  perpendicular 
to  experimental  interstage  plate  -  container  axis  -  degrees 

Radius  of  the  spherical  jet  boundary  or  floor  shock  approxi¬ 
mations  -  radii 

Angle  between  a  line  of  experimental  pressure  taps  and  a  line 
from  the  nozzle  centerline  perpendicular  to  a  splitter  fence, 
measured  positive  counter-clockwise  looking  toward  nozzle  - 
degrees 

Angle  between  sides  (a)  of  first  stagnation  control  volume 
(V)  -  see  Figure  (8) 

Angle  between  jet  axis  and  perpendicular  to  the  interstage  floor 
shock  -  degrees 

Angle  between  a  line  of  experimental  pressure  taps  and  the  posi¬ 
tive  right-hand  vector  direction  of  nozzle  tilts  ( {$  ),  measured 
clockwise  looking  at  nozzle  -  degrees 
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in  the  ambient  air  or  in  the  interstage  cavity 
Viilue  assumed  at  the  start  of  an  iteration 
at  the  separation  rocket  jet  boundary 

at  the  separation  rocket  jet  centerline,  or  referring  to  mass 
flow  from  nozzle 

value  computed  at  the  end  of  an  iteration 
at  the  exit  of  the  first  stagnation  layer 
at  a  point  on  the  interstage  floor 
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_ Explanation _ 

at  a  point  in  the  separation  rocket  free  jet 

at  a  point  in  the  first  stagnation  layer 

at  the  separation  rocket  nozzle  exit 

at  a  point  on  the  experimental  nozzle  plate 

between  the  jet  and  experimental  plate  -  container  centerline 
at  the  level  of  the  plate 

out  of  a  volume 

at  a  point  on  the  'experimental  interstage  floor  plate 
at  a  point  on  an  experimental  pitot  -  tube  rake 
in  the  jet  radial  direction 
Total  or  stagnation  value  in  a  flow 

in  the  jet  axial  direction,  or  referring  to  flux  of  axial  momentum 
from  nozzle 

at  the  nozzle  throat 
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INTRODUCTION 


Nature  of  the  Recirculation  Problem 

Recirculation  flow,  causing  gas  dynamic  interactions  between  separating  vehicle 
stages,  is  common  when  staging  is  at  high  altitude  and  the  separating  stage  is 
powered  by  a  rocket  motor  which  is  ignited  during  the  staging  sequence.  Under 
these  conditions  the  jet  wake  from  the  rocket  motor,  expanding  to  the  low  ambient 
pressure,  may  fill  the  space  between  the  stages  and  stagnate  directly  on  the  floor 
of  the  interstage  volume.  When  this  occurs  (Figure  l)  a  strong  peak  in  floor 
pressure  is  set  up  where  the  jet  flow  is  normal  to  the  surface,  with  pressure  de¬ 
clining  rapidly  away  from  the  peak  in  all  directions.  These  floor  pressure 
gradients  cause  the  stagnating  flow  to  re -accelerate  to  supersonic  speeds  laterally 
along  the  floor  surface  away  from  the  peak.  A  second  stagnation  may  take  place  with 
supersonic  flow  recirculating  up  toward  the  separating  stage,  if  there  is  any  inter¬ 
stage  hardware  blocking  lateral  escape  of  gas  (Figure  la),  or  if  there  is  an  adjacent 
jet  stagnation  flow  field  opposing  the  lateral  flow  of  gas  from  under  the  separating 
stage  (Figure  lb).  Recirculation  is  the  existence  of  any  supersonic  reverse  flow 
directed  back  onto  the  separating  stage,  which  arises  from  the  deflection  of  exhaust 
gas  from  the  stage’s  own  rocket  motor. 

Although  there  may  be  substantial  gas  dynamic  interactions  between  the  stages,  recir¬ 
culation  does  not  occur  during  fire-in-the-hole  separation  (initially  closed  intar- 
stage),  such  as  Polaris  first  separation.  This  is  because  a  high  pressure  is  rapidly 
achieved  in  the  interstage  volume  which  confines  the  separation  rocket  jet  wake  and 
suppresses  supersonic  stagnations  and  re-expansions.  Supersonic  recirculation  flow 
gradually  becomes  important  as  more  and  more  venting  is  allowed  in  the  interstage 
design,  provided  the  separation  occurs  at  high  altitude,  so  that  ambient  pressure 
does  not  confine  the  flow  in  the  interstage.  Because  the  interstage  pressure  must 
be  high  to  suppress  recirculation,  fire-in-the-hole  separations  are  characterized 
by  higher  staging  forces  than  separations  in  which  recirculation  predominates. 
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However,  when  recirculation  is  present,  in  a  vented  design,  staging  forces  may  be 
much  higher  than  would  be  expected  assuming  complete  mixing  within  the  interstage  for 
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that  design#  This  is  because  the  momentum  of  the  reverse  flow  jet  is  available  to 
act  on  the  separating  stage,  in  addition  to  the  interstage  cavity  pressure,  and  for 
supersonic  reverse  flow,  this  momentum  may  be  quite  large* 

Importance  of  Recirculation  to  Vehicle  Designs 

The  presence  of  supersonic  recirculation  flow  within  an  interstage  usually  poses 
serious  vehicle  design  and  performance  problems.  Often  the  principle  problem  is 
the  large  force  produced  on  the  separating  3tage  base  by  the  stagnation  of  the  re¬ 
verse  flow  jet.  The  axial  component  of  this  force  imparts  an  unexpected  velocity 
to  the  separating  stage,  changing  its  trajectory.  If  the  separating  3tage  is  unguxded, 
no  in-flight  correction  can  be  made  for  this  velocity  change.  Even  in  guided  stages, 
the  velocity  perturbation  often  occurs  in  too  short  a  time  to  be  sensed  accurately  by 
the  guidance  system.  This  often  results  in  vehicle  accuracy  problems.  For  example, 
the  recirculation  occurring  during  separation  of  the  unguided  A3  re-entry  system  shifts 
the  point  of  impact  down  range  by  approximately  six  miles. 

Because  the  recirculation  flow  is  supersonic,  the  pressures  it  produces  on  both  stages 
may  be  unsymmetrical,  imparting  angular  rates  to  the  stages.  This  occurs  if  the  inter¬ 
stage  hardware  is  unsyiranetrical,  or  if  adjacent  separating  stages  have  different  ig¬ 
nition  times  or  different  thrust  levels.  Angular  rates  during  staging  may  cause  clearance 
problems,  and  if  the  stages  are  unguided,  their  trajectories  will  be  affected.  For 
example,  the  position,  orientation,  and  radar  cross  section  of  the  Polaris  A3  re-entry 
system  components  have  been  drastically  affected  by  recirculation-induced  angular  rates, 
and  during  re-entry,  these  rates  cause  re-entry  body  limit-cycle  wobbling  oscillations 
which  increase  radar  scattering  by  the  wake. 

Finally,  important  design  problems  are  often  caused  by  unexpectedly  high  structural 
loads  and  heating  produced  in  the  interstage  by  the  recirculation  flow  field.  For 
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example,  the  high  loads  experienced  on  the  Polaris  A3  forward  heat  shield,  during 
re-entry  body  separation,  iiava  required  about  50  lb.  of  extra  structure  in  this 
area.  Even  so,  the  smaller  components  (firing  units,  cables,  snubbers,  etc.)  are 
regularly  torn  off  the  missile  by  the  recirculation  flow  field  and  scattered  about 
the  sky.  This  has  greatly  complicated  the  problem  of  radar  data  analysis  from  test 
flights. 

Expected  Future  Occurrence  of  Recirculation 

Recirculation  problems  during  re-entry  system  separation  on  the  Polaris  A3  will  con¬ 
tinue  to  be  important  as  the  system  undergoes  modifications.  Even  minor  changes,  such 
as  different  instrumentation  layouts  in  test  flights,  have  been  observed  to  change  the 
recirculation  effects  substantially.  Large  changes,  such  as  lnvisioned  for  the  Mark 
XU  adaptation,  will  probably  introduce  completely  new  recirculation  problems.  Super¬ 
sonic  flow  in  the  interstage  which  has  been  shown  to  occur  during  B3  first  stage 
separation  (Reference  a)  must  also  be  solved. 

Recirculation  effects  should  be  anticipated  in  any  highly  vented  interstage  design, 
such  as  envisioned  in  various  large  solid  booster  concepts.  For  example,  interstage 
recirculation  gas  dynamics  produce  important  problems  in  the  highly-vented  Saturn 
interstage  (Reference  b).  Recirculation  forces,  and  debris  carried  by  the  recirculation 
flow  field,  may  be  important  during  lunar  landing  and  take-off  (Reference  c),  especially 
if  there  is  substantially  irregularity  in  the  landing  surface. 

Available  Methods  and  Data  for  Recirculation  Problems 

The  flow  in  the  expanding  jet  from  the  separation  rocket  nozzle  may  be  calculated  by 
the  LMSC  method-of-characteristics  computer  program  of  Reference  (d).  A  Newtonian 
approximation  to  the  interstage  floor  pressures,  produced  by  the  direct  stagnation  of 
this  jet,  may  be  calculated  with  the  jet  impingement  computer  program  of  Reference  (e). 
However,  these  pressures  may  be  substantially  in  error  if  a  recirculation  flow  field 
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exists  between  the  free- jet  flow  and  the  interstage  floor  surface.  When  the  inter¬ 
stage  is  initially  closed  ( fire-in-the-hole ) ,  so  that  recirculation  is  suppressed, 
a  computer  program  assuming  complete  mixing.  Reference  (f),  may  be  used  to  calculate 
the  axial  force  during  separation.  However,  this  method  gives  no  information  on 
local  flow  conditions  within  the  interstage.  Also,  it  shows  poor  agreement  with 
experiment  when  the  separation  plane  is  near  the  interstage  floor,  or  where  venting 
is  initially  present  near  the  floor.  When  recirculation  flow  is  responsible  for  the 
staging  forces,  there  is  no  analytical  method  for  predicting  magnitudes  or  determining 
local  flow  conditions  over  interstage  hardware. 

Large  amounts  of  experimental  data  on  Polaris  A3  re-entry  body  separation  forces  and 
interstage  pressures  have  been  accumulated  at  LMSC  (References  (g),  (h),  and  (i)). 

These  data  are  for  reasonably  exact  scale  models  of  different  re-entry  system  designs, 
and  consequently  the  geometries  are  very  complicated.  Further,  pressure  data  are  very 
limited,  with  no  pitot  tube  or  bther  surveys  of  recirculation  flow  conditions.  Con¬ 
siderable  interstage  static  pressure  and  vent  flow  pitot  pressure  data  have  been  taken 
on  the  separation  of  a  small  body  from  a  highly  vented  container  (Reference  j),  in 
which  the  vents  are  longitudinal  and  not  axi symmetric,  and  the  interstage  floor  may 
have  a  shape  which  varies  *Jith  time  in  a  complicated  way.  Reference  (b)  contains 
interstage  floor  pressure  data  for  designs  with  similar  vents  to  those  of  Reference  (j), 
and  various  (fixed)  floor  shapes.  There  are  no  flow  property  surveys  in  Reference  (b). 
Reference  (k)  presents  pressure  data  on  a  flat  plate  due  to  the  impingement  of  various 
jets  at  low  ambient  pressures,  and  gives  Schlieren  photographs  showing  jet  boundaries 
and  stagnation  shock  shapes. 

Objectives  of  the  Independent  Development  Study 

Although  attempts  have  been  made  to  develop  analytical  methods  for  predicting  recircu¬ 
lation  effects,  using  the  available  experimental  data  (especially  that  of  Reference  i), 
the  complications  of  the  configurations  tested  have  proved  insurmountable.  To  develop 
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an  analytical  capability,  consideration  of  simple  configurations  ia  necessary,  with 
subsequent  generalisations  to  more  complicated  designs*  However,  no  experimental 
data  exists  on  simple  interstages,  with  which  to  verify  a  theoretical  analysis  (except 
for  Reference  (k)  in  which  there  ia  no  interstage  at  all)*  The  importance  of  recir¬ 
culation  estimates  in  the  design  of  advanced  Polaris  re-entry  systems,  and  the  lack 
of  fundamental  theoretical  knowledge  about  recirculation,  led  to  the  present  study  on 
196U  Augmented  Independent  Development  Funds,  to  develop  a  basic  understanding  of  the 
recirculation  process* 

The  objectives  of  the  Qas  Dynamics  Interactions  study  were: 

1)  To  develop  an  analytical  method  for  computing  recirculation  staging 
forces  for  simple  designs  in  which  the  interstage  is  symmetric  about 
the  separation  rocket  nozzle  axis* 

2)  To  generalize  the  axisymmetric  method  to  various  idealized  types  of 
asymmetries  found  in  actual  vehicle  designs* 

3)  To  accumulate  parametric  experimental  data  on  recirculation  flow  field 
properties  and  induced  staging  forces,  for  idealized  axisymmetric  and 
ncu~ax± symmetric  designs,  and  use  these  data  to  assist  in  the  development 
of  the  theory  and  provide  information  on  trends  which  could  be  used  in 
design  estimates* 

Because  of  the  emphasis  on  idealized  configuration  in  both  the  theoretical  and  experi¬ 
mental  investigations,  this  work  was  supported  by  Independent  Development  Funds.  By 
divorcing  the  study  from  any  specific  vehicle  design,  it  was  hoped  the  results  would 
apply  to  preliminary  design  estimates  for  many  different  types  of  vehicles* 
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PROJECT  PLAN 


Original  Plan  and  Funding 

Theoretical  analysis  was  planned  for  a  series  of  progressively  more  complicated  axi- 
symmetric  interstage  designs,  beginning  with  the  first  stagnation  of  the  separation 
rocket  jet  on  ’.fates  representing  interstage  floors  of  various  shapes,  then  proceeding 
to  add  tho  sides  and  top  of  the  interstage,  and  to  consider  various  types  of  venting 
(see  figure  2).  Generalization  of  the  axisymmetric  method  was  to  be  attempted  for 
i  .olination  of  the  separating  stage  from  the  interstage  axis,  and  for  the  presence 
of  adjacent  separating  stages  assumed  to  be  identical  to  the  stage  under  consideration# 
Experimental  data  was  to  be  taken  on  each  of  the  idealized  configurations  considered 
in  the  analysis,  so  as  to  support  and  verify  each  step  of  the  theory.  Sufficient 
pressure  distributions,  flow  field  property  surveys,  and  photographs  of  flow  field 
structures  were  tc  be  taken  to  permit  thorough  checkout  of  the  assumptions  made  in 
the  theoretical  work. 

The  plan  called  for  testing  to  begin  as  soon  as  possible  after  approval  of  the  project. 
To  accumulate  large  amounts  of  data  quickly  and  at  low  cost,  cold-flow  testing  was 
planned  for  LMSC's  Tenney  Altitude  Chamber  (8U— 33 ) ®  Design  and  fabrication  of  the 
simple  axisymmetric  model  configurations  was  scheduled  to  be  complete  by  mid-August, 
with  the  first  data  to  be  available  for  analysis  by  the  end  of  August.  Testing  of 
more  complicated  models  was  to  continue  to  the  end  of  the  project. 

The  development  of  a  theoretical  method  for  computing  recirculation  forces  in  the 
axisymmetric  design  was  planned  to  begin  before  the  testing,  and  to  be  complete  by 
the  end  of  October,  allowing  two  months  for  the  use  of  test  data  in  this  effort. 
Generalization  of  the  axisymmetric  method  to  the  various  asymmetries  was  to  begin  in 
November  and  to  continue  concurrent  with  testing  of  the  corresponding  configurations. 

The  total  budget  for  the  project  was  $95,000.  Of  this  amount,  $$%  was  allocated  to 
engineering  work  on  the  theory,  data  correlations,  model  design,  and  test  conduct* 
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Approximately  3 0^  oft  the  budget  was  to  be  used  for  facility  charges  and  data  re¬ 
duction  computer  programming.  1$%  of  the  budget  was  for  model  fabrication. 

Fulfillment  of  Project  Plan 

The  objectives  of  the  experimental  phase  were  met  in  full.  Testing  was  completed  on 
all  the  configurations  planned,  with  several  configurations  and  conditions  tested 
which  were  not  originally  scheduled.  Adequate  pressure  surveys,  loads  and  photographic 
data  were  obtained  to  permit  theoretical  analysis  of  all  the  configurations.  However, 
the  start  of  testing  was  delayed  to  the  end  of  September  because  of  lack  of  personnel 
available  at  the  beginning  of  the  project,  problems  with  facility  scheduling,  and 
problems  with  the  photographic  setup.  The  first  data  were  not  available  for  analysis 
until  the  beginning  of  November.  The  test  schedule  was  very  much  compressed  to  finish 
the  planned  testing  within  allocated  funding.  Consequently,  much  data  was  accumulated 
with  very  little  concurrent  analysis  of  the  results. 

The  theoretical  analysis  of  the  first  stagnation  of  the  separation  rocket  jet  on  the 
interstage  floor  was  completed  by  the  end  of  August,  and  results  from  an  RPC  1000  com¬ 
puter  program  were  obtained  by  the  end  of  September.  Calculation  of  further  flow  re¬ 
accelerations  and  stagnations  within  the  axisymne trie  interstage  had  been  studied,  and 
an  RPC  1000  program  had  been  drawn  up  by  the  end  of  October.  With  the  availability  of 
experimental  data  in  November,  it  was  found  that  the  First  Stagnation  program  failed 
to  agree  with  experiment  for  many  cases  of  interest.  The  remaining  analysis  time  was 
spent  in  studying  the  data  to  determine  why  the  First  Stagnation  theory  was  failing, 
how  it  could  be  corrected,  and  what  implications  could  be  drawn  concerning  the  inter¬ 
stage  calculation  as  a  whole.  No  progress  was  made  on  the  calculation  of  the  various 
non-axi symmetric  test  configurations. 
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DESCRIPTION  0?  EXPERIMENTAL  WORK 


Model  Description 

The  test  was  conducted  at  the  IMSC  Tenney  Altitude  Facility  during  the  period  September- 
October  1 A  schematic  of  the  test  set  up  is  presented  in  Figure  (3)*  High  pressure 
nitrogen  gas  fran  12  standard  ccnmercial  bottles,  manifolded  in  parallel,  flowed  through 
a  valve  and  out  through  a  nozzle  onto  a  plate  or  into  a  vented  container  within  the 
evacuated  test  chamber.  Thrust  produced  by  the  nozzle  and  the  force  acting  on  the  plate 
were  measured  using  load  cells*  Pressures  measured  included  ambient  (test  chamber) 
pressure,  nozzle  chamber  pressure,  static  pressures  on  the  plate  and  on  the  walls  of 
the  containers,  and  pitot  pressures  in  the  jet  wake  and  in  the  vent  areas.  High  speed 
Schlieren  motion  pictures  were  taken  during  most  runs.  The  axial  distance  between  the 
nozzle  and  the  plate  and  the  container  vent  area  were  varied  from  run  to  run. 

To  facilitate  completing  the  test  schedule  within  the  allotted  time,  the  test  set-up 
was  planned  for  maximum  convenience  of  operation.  Nitrogen  pressure  regulation,  valve 
control,  and  camera  frame  rate  and  on-off  control  were  accomplished  outside  the  test 
chamber.  The  separation  distance  between  the  nozzle  and  all  the  various  plates  or 
containers  was  remotely  controlled  using  an  electric  motor  to  drive  a  worm  and  wheel 
gear  arrangement. 

Two  different  nozzles  were  used  during  the  test  program:  a  sonic  nozzle  and  a  super¬ 
sonic  nozzle  having  an  expansion  ratio  of  £.5  to  1.  These  nozzles  are  shown  in  Figure 

(U). 

As  shown  in  Figure  (£),  interstage  floor  configurations  tested  included  a  flat  plate, 
and  a  conical  plate.  Certain  model  configurations  included  a  flat  plate  concentric 
with  the  nozzle  axis  and  mounted  in  the  plane  of  the  nozzle  exit  (nozzle  plate).  Several 
runs  were  made  with  the  nozzle  tilted  with  respect  to  the  plate-container  axis.  A  cir¬ 
cular  fence  two  inches  high  and  four  inches  in  radius  mounted  concentrically  about  the 
noizle  axis,  was  used  with  the  flat  plate  to  simulate  a  vented  container  into  wM  tbs 
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Jet  exhausted*  The  fence  was  moved  axially  between  the  plate  and  the  nozzle  to 
furnish  different  vent  areas.  The  circular  fence  waa  also  used  with  the  nozzle 
plate  to  simulate  a  re-entry  body  flare.  The  effeot  of  adjaoent  jets  was  simu¬ 
lated  using  a  flat  plate  and  splitter  fences  (see  Figure  5).  The  preaenoe  of 
two,  throe,  and  six  jets  was  thus  simulated. 

The  parameters  that  were  varied  from  run  to  run  (see  Figures  6  and  7)  included 
the  separation  distance  between  the  plate  and  the  nozzle  exit  plane  (x),  the  dis¬ 
tance  from  the  plate  to  the  bottom  of  the  circular  fence  (h),  and  the  radial  dis¬ 
tance  between  the  nozzle  axis  and  the  splitter  fences  (f).  The  plate  and  the 
fences  were  attached  so  that  they  had  the  same  movement  relative  to  the  nozzle  exit 
plane. 

On-off  control  of  the  nitrogen  flow  was  through  a  solenoid  controlled,  hydraulically 
actuated  quick  opening  Jamesbury  ball  valve.  Constant  nitrogen  pressure  during  a 
run  was  maintained  with  a  dome  regulator. 

Test  Conditions 

Before  every  test  was  run  the  Tenney  Chamber  was  evacuated  to  a  pressure  altitude  of 
approximately  117,000  feet.  After  a  test  run  (each  run  lasted  one  second)  the  altitude 
had  dropped  to  about  105,000  feet,  An  increase  in  test  chamber  pressure  of  about  65 
percent.  The  nozzle  chamber  pressure  averaged  about  11*50  psi  during  a  run.  This 
average  chamber  pressure  varied  from  run  to  run,  being  dependent  on  the  bottle  pressure. 
Approximate  extremes  of  avorage  chamber  pressure  were  1350  psi  and  1510  psi.  Within  a 
test  run  the  cliamber  pressure  reached  operating  level  within  30  milliseconds  after  the 
valve  was  actuated  and  would  drop  approximately  50  psi  per  second  in  most  cases.  The 
epeed  at  which  nozzle  pressure  decreased  (or  increased)  depended  on  the  pressure  in  the 
bottles. 


As  shown  in  Figure  (6),  separation  distance  (x)  varied  from  10.7  r„  (nozzle  exit  radii) 


to  h9«$  when  the  sonic  nozzle  was  used,  and  from  to  19*25  when  the 

eupersonic  nozzle  was  used.  Cylindrical  fence  height  was  2n  (21*7  i'Vi  sonic,  9*07 
supersonic)*  The  distance  from  the  plate  to  the  bottom  of  the  circular  fence  (h) 
varied  from  0  to  21*7»hSonic  and  from  0  to  U*53^  supersonic*  Maximum  lower  vent 
area  formed  by  this  separation  between  plate  and  fence  was  7,8UO  A*  (nozzle  throat 
area)  for  the  single-stage  configurations,  and  with  the  adjacent  jet  (splitter  fence) 
models  the  maximum  area  was  I960  Aft*  Only  the  sonic  nozzle  was  used  with  the  splitter 
fences* 

Instrumentation 

Instrumentation  used  in  the  test  program  included  strain  gage  pressure  transdrcers  and 
load  cells,  a  position  transducer,  thermo-couples,  bourdon  tube  pressure  gages  for 
measuring  nitrogen  supply  pressure,  and  a  bourdon  tube  vacuum  gage  for  measuring  test 
chamber  pressure. 

Thrust  produced  by  the  nozzle,  and  the  force  produced  by  the  jet  impinging  on  the  plate 
or  container,  were  measured  using  single-bridge  strain  gage  load  cells* 

Strain  gage  pressure  transducers,  of  both  absolute  and  differential  types,  were  used 
to  measure  test  chamber  pressures  and  model  pressures,  including  pitot  pressures  in 
the  jet  wake  and  in  vent  areas  and  static  pressures  on  the  plates  and  fences.  Typical 
locations  for  these  measurements  are  shown  in  Figure  (6),  A  position  transducer,  the 
output  of  which  was  on  digital  display,  was  used  to  set  the  separation  distance.  Two 
thermocouples  were  used  to  measure  temperatures  in  the  shock  layer  on  the  plate.  In 
addition  to  force  and  pressure  instrumentation,  extensive  use  was  made  of  photographic 
instrumentation  consisting  of  a  modified  schlieren  system  and  a  high  speed  motion 
picture  camera.  As  shown  in  Figure  (3),  light  from  a  point  source  was  passed  through 
the  test  model  and  returned  to  the  camera  lens  along  the  same  path  by  a  12"  diameter 
aBpheric  mirror* 
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All  load  cells  and  pressure  transducers  were  physically  calibrated  before  their  in¬ 
stallation  in  the  model*  In  this  physical  calibration  the  output  of  the  transducer 
was  read  on  a  strain  indicator*  Also,  with  no  load  on  the  transducer,  known  re¬ 
sistances  were  shunted  into  the  circuit,  producing  changes  in  the  strain  indicator 
reading.  Thus,  stops  were  obtained  relating  apparent  load  on  the  transducer  with 
resistance  in  the  circuit.  These  steps  were  shunted  into  the  circuit  to  make  the 
electrical  transducer  calibratior^i  before  each  test  run* 

Data  Acquisition 

Analog  signals  from  the  strain  gage  load  cells  and  pressure  transducers  were  digitized 
and  recorded  by  the  Data  Acquisition  System  (DAS)  of  0/8l;-33»  A  digitizer  accepted 
analog  input  voltages  from  the  transducers  and  load  cells  and  derived  a  number  pro¬ 
portional  to  the  amplitude  of  each  voltage  at  the  time  of  sampling.  This  number 
(millivolts)  was  then  written  on  magnetic  tape  In  BCD  format*  Time  identification  was 
recorded  with  the  test  data*  After  each  series  of  runs,  the  data  was  dumped  onto  cards 
and  also  listed  by  an  IBM  I4O7  accounting  machine* 

Before  each  test  run,  with  the  test  chamber  at  the  proper  simulated  altitude,  no-load 
transducer  outputs  were  recorded,  then  an  electrical  calibration  made*  In  this  electri¬ 
cal  calibration  a  known  resistance  (the  same  as  used  in  the  pre-test  physical  calibration 
was  shunted  into  each  transducer  circuit,  producing  a  change  in  millivolt  counts  in  the 
DAS*  This  permitted  deriving  a  relationship  between  load  (force  or  pressure)  on  the 
transducer  and  electrical  output  in  millivolt  counts* 

Test  Summary 

The  test  program  required  approximately  ten  weeks  to  complete,  including  three  weeks  for 
calibration  and  installation  of  transducers  and  set-up  of  the  model  and  the  nitrogen 
supply*  During  the  seven  weeks  of  testing,  2 $9  runs  were  made,  167  with  the  sonic  nozzle 
and  92  with  the  supersonic  nozzle* 


A  complete  summary  of  the  experimental  configurations  may  be  found  in  Tables  I,  n 
and  HI,  Table  I  summarizes  the  sonic  nozzle  configurations  tested,  while  Table  II 
presents  the  supersonic  nozzle  configurations.  Table  HI  summarizes  the  adjacent 
jet  Simula ting  fence  configurations,  all  of  which  were  tested  with  the  eonio  nozzle* 
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EXPERIMENTAL  RESULTS 


Data  Reduction 

To  facilitate  the  reduction  of  the  data,  obtained  as  counts  (see  previous  section), 
a  709k  Fortran  data  reduction  computer  program  was  written  (Reference  r).  This  pro¬ 
gram  was  designed  to  quantify  and  list  the  experimental  data.  The  program  converts 
as-received  counts  into  data  expressed  in  engineering  units.  This  data  is  listed  as 
a  function  of  time,  as  is  normalized  data.  Data  is, normalized  to  a  nozzle  source 
pressure  of  1000  psia  by  dividing  the  data  by  the  instantaneous  average  nozzle 
pressure  and  multiplying  the  result  by  1000.  In  addition,  the  program  permits  calcu¬ 
lation  of  average  values  of  the  data  over  a  specified  time  interval,  as  well  as 
generation  of  plots  of  any  or  all  data  versus  time  for  the  same  time  interval. 

IXie  to  the  fact  that  some  of  the  pressure  transducers  measured  differential  pressure, 
while  others  measured  absolute,  additional  operations  were  performed  on  the  reduced 
data  to  give  all  pressure  measurements  in  absolute  units.  To  obtain  the  absolute 
loading  on  both  the  flat  plate  and  nozzle  plate  -  nozzle  combination,  it  was  necessary 
to  convert  the  differential  force  measured  experimentally  by  adding  the  ambient 
pressure  loads  acting  on  the  test-cell  sides  of  the  plates.  The  total  normalized 
force  in  each  case  was  then  divided  by  the  average  normalized  nozzle  thrust  taken 
from  all  the  flat  plate  test  runs. 

The  average  nozzle  thrust  was  used  to  determine  the  force/thrust  ratio  rather  than  the 
thrust  for  the  test  run  in  question,  for  two  reasons.  First,  for  runs  which  included 
the  nozzle  plate,  the  measured  force  on  the  nozzle  plate  included  the  pressure  loading 
on  the  plate  due  to  recirculation  as  well  as  the  nozzle  thrust.  Since  there  was  no  way 
to  accurately  separate  the  two  contributions  to  the  force,  the  nozzle  thrust  for  those 
runs  was  unknown.  Second,  a  plot  of  nozzle  total  pressure  versus  measured  -thrust  did 
not  reveal  the  expected  linear  relationships.  Since  for  a  given  nozzle,  the  thrust 
should  be  directly  proportional  to  the  nozzle  total  pressure,  other,  unknown  effects 
made  their  presence  felt.  These  effects  may  have  included,  for  example,  errors  in 


force  measurement  or  nozzle  flow  problems,  such  as  the  effective  half  angle  or  exit 
Mach  number  differing  from  design  conditions.  The  value  of  average  thrust  was  35.U95 
pounds  for  the  sonic  nozzle,  and  hl*hS7  pounds  for  the  supersonic  nozzle*  Both  of 
these  figures  are  normalized  to  a  nozzle  total  pressure  of  1000  paia. 

Description  of  Experimental  Data 

For  purposes  of  presentation,  the  data  was  divided  into  three  classifications  as 
follows: 

1*  Sonic  Nozzle  Configurations 
2,  Supersonic  Nozzle  Configurations 
3*  Adjacent  Jet  Simulating  Fence  Configurations 
A  summary  of  the  configurations  tested  can  be  found  in  Tables  I/,  II,  and  HI* 

Plots  of  the  experimental  test  data  for  the  first  two  classifications,  sonic  and  super¬ 
sonic  nozzle  configurations,  may  be  found  in  Appendix  II.  The  nozzle  total  pressure 
to  ambient  pressure  ratio  for  the  data  presented  in  the  appendix  is,  nominally,  12,000. 

In  Appendix  II  plots  of  plate  pressure  distributions,  nozzle  plate  pressure  distributions, 
fence  pressure  distributions,  plate  force/thrust  ratios  vs  separation  distance,  nozzle 
plate  force/thrust  ratios  vs  separation  distance,  and  rake  pressures  vs  separation  dis¬ 
tance,  are  presented.  The  photographic  data  on  flow  field  structures  is  not  presented 
due  to  lack  of  time  available  to  scale  and  map  the  photographs  to  compensate  for  angle- 
of-view  distortions  (see  Comparison  of  First  Stagnation  Theory  to  Experiment,  for  a 
discussion  of  the  analysis  required). 

Due  to  time  limitations,  the  data  from  the  third  classification,  the  adjacent  jet  simu¬ 
lating  fence  configurations,  is  not  included  in  Appendix  II.  However,  these  data  were 
reduced  and  rough  plots  made  during  the  data  analysis  effort,  so  that  the  data  are 
available  for  use  in  theoretical  work  and  design  estimation. 

Data  Validity 


To  assess  the  validity  of  the  data,  a  comparison  of  the  flat  plate  pressure  distributions 
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from  Appendix  II  with  the  experimental  data  of  Reference  (k),  and  with  the  theoretical 
predictions  of  the  program  of  Reference  (e)  has  been  made  in  Figure  (lii)# 

For  purposes  of  comparison,  it  was  necessary  to  cross  plot  tne  experimental  data  of 
Stitt  (Reference  k)  to  obtain  data  at  the  same  separation  distances  as  were  tested* 
Further,  it  should  be  noted  that  Stitt’s  data  was  obtained  at  a  pressure  ratio  of  1*52 
x  10^  while  the  present  data  holds  for  a  ratio  of  1*2  x  10^*  This  difference  affects 
the  tails  of  the  distribution,  since  the  pressure  ratio  determines  the  boundary  of  the 
jet*  The  core  region  does  not  appear  to  be  greatly  affected  by  the  ratio,  however* 

It  can  be  seen  that  the  experimental  data  from  the  test  and  from  Stitt  are  in  good 
agreement. 

>*  •  t 

The  Newtonian  jet  impingement  theory  of  Reference  (e),  compared  in  Figure  (lij),  assumes 
the  plate  shock  lies  directly  on  the  surface,  with  no  standoff  distance.  This  would 
be  expected  to  give  lower  pressures  on  a  flat  plate  than  actually  occur  because  the 
Mach  number  in  the  jet  is  higher  at  the  plate  than  at  any  shock  standoff  distance 
above  the  plate.  This  is  borne  out  in  Figure  (lli),  where  jet  impingement  data  is  lower 
at  the  axis  than  either  experimental  data  source.  The  magnitude  of  this  difference 
is  consistent  with  the  standoff  distance  observed  in  the  photographs  of  the  test 
configurations* 


DESCRIPTION  OF  THEORETICAL  WORK 
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Introduction 

As  described  under  Project  Plan,  the  theoretical  studies  were  concentrated  on  a 
series  of  progressively  more  complicated  axially  symmetric  interstage  designs, 
beginning  with  the  problem  of  a  jet  stagnating  on  a  plate,  then  considering 
stagnations  against  the  sides  and  top  of  the  interstage*  This  organization 
stemmed  from  preliminary  concepts  concerning  interstage  recirculation  gas  dynamics 
drawn  from  jet  impingement  studies  (e*g.  Reference  c)  and  the  data  of  References 
(i)  and  (k).  It  appeared  that  the  flow  in  the  interstage,  for  recirculation  cases, 
was  a  step-wise  process  of  stagnations  and  re-accelerations,  with  each  stagnation 
process  separated  from  the  preceding  one  by  a  supersonic  region  (see  Figure  la)* 

Thus  it  was  thought  that  each  area  of  stagnation  could  he  solved  independently 
from  the  other  areas,  given  only  the  interstage  pressure,  the  interstage  geometry, 
and  the  supersonic  flow  conditions  upstream  of  the  stagnation.  By  solving  the 
stagnation  of  the  jet  on  the  plate  to  obtain  the  flow  conditions  leaving  the  stag¬ 
nation  layer,  it  was  thought  that  the  initial  conditions  would  then  be  given  for 
calculation  of  the  acceleration  of  gas  along  the  plate  and  stagnation  against  the 
interstage  sides.  Likewise,  solution  of  the  side  stagnation  was  anticipated  to 
set  up  the  independent  problem  of  the  stagnation  against  the  top  of  the  interstage* 

A  scheme  for  solving  the  axisymmetric  interstage  was  envisioned  in  which  interstage 
cavity  pressure  (Pc*  of  Figure  1)  would  be  first  assumed,  then  the  chain  of  stagnations 
and  re-accelerations  areound  the  interstage  would  be  calculated,  and  from  the  degree 
to  which  the  mass  flow  out  of  the  interstage  found  for  the  various  exits  was  different 
from  the  mass  flow  in  from  the  nozzle,  a  new  interstage  pressure  would  be  assumed.  The 
process  would  be  repeated  until  the  interstage  was  in  balance.  Thus,  the  theoretical 
analysis  of  the  first  stagnation  of  the  separation  rocket  jet  on  a  plate,  given  the 
plate  shape  and  the  surrounding  ambient  (or  cavity)  pressure,  was  considered  as  an 
independent  problem,  with  application  of  the  results  to  the  re-acceleration  and 
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side  stagnation  a  separate  study# 

As  has  been  indicated  in  Project  Plan,  the  analysis  of  the  first  stagnation  was 
completed,  but  based  on  assumptions  which  make  the  analysis  inadequate  in  the  high- 
pressure  confined  environment  of  the  interstage*  However,  this  method  does  show 
reasonable  agreement  with  experiment  in  very  low  ambient  pressure  cases,  and  might 
be  modifiable  for  high  pressure  cases.  Consequently,  a  discussion  of  this  work 
is  given  below.  The  analysis  of  the  side  stagnation  was  set  up,  but  without  good  flow 
conditions  predicted  from  the  first  stagnation  within  the  interstage,  the  adequacy 
of  the  side  stagnation  theory  could  not  be  investigated.  Theoretical  work  was  done 
on  first  stagnation  for  the  nozzle  tilted  with  respect  to  the  plate,  but  this  study 
could  not  be  carried  far  enough  to  produce  results. 

First  Stagnation  Theory 

Two  methods  were  investigated  for  calculating  the  normal  stagnation  of  a  jet  on  a 
concentric  ari  symmetric  plate.  A  stream  tube  continuity  approach  along  the  lines 
suggested  by  Reference  (m)  was  attempted,  which  would  give  the  detailed  flow  field 
within  the  stagnation  layer.  Persistent  mathematical  difficulties  were  encountered 
and  a  working  method  was  not  completed.  A  second  approach  was  to  make  a  control- 
volume  analysis  of  the  layer  (volume  7  of  Figure  8),  to  determine  the  total  load  on 
the  plate,  the  height  of  the  layer  at  the  jet  boundary,  and  the  average  flow  con¬ 
ditions  crossing  the  exit  of  the  layer  (d  of  Figure  8).  This  method  was  considered 
adequate  to  solve  the  staging  force  problem  and  pro.ide  input  flow  conditions  to  the 
side  stagnation  calculation#  A  working  method  and  an  RPC  1*000  computer  program  were 
developed  using  the  control -volume  analysis. 

A  full  discussion  of  the  control-volume  method  is  given  in  Appendix  I.  The  basic 
equations  are  conservation  of  mass,  axial  momentum,  and  radial  momentum  for  volume 
7  of  Figure  (8).  These  equations  introduce  (?)  unknown  quantities  and  thus  (6) 
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additional  relationships  are  required.  The  basic  assumptions  of  the  method  are  that 
the  jet  boundary  is  given  by  the  circular  arc  method  of  Reference  (n),  the  stagnation 
shock  is  a  spherical  cap,  and  only  isentropic  jet  flow  as  given  by  the  program  of 
Reference  (d)  crosses  the  stagnation  shock  (i.e.  the  thickness  and  strength  of  the 
jet  boundary  region  is  negligible)*  The  shock  stand  off  distance  on  the  Jet  center¬ 
line  is  taken  from  Reference  (o).  These  assumptions  permit  relations  to  be  drawn 
between  the  unknown  layer  exit  height  (  of  Figure  8)  and  exit  radius  (  ry  ), 

exit  average  flow  angle  (  ),  exit  total  pressure  (fu,  )  radial  momentum  flux 

entering  the  layer,  and  layer  cross-sectional  area.  The  sixth  relation  is  the 
assumption  that  average  layer  static  pressure  equals  average  plate  pressure.  With 
these  relations,  the  three  basic  equations  are  solved  by  a  doable-iteration  scheme. 
Topical  results  of  the  analysis  are  shown  in  Figures  9-11* 
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COMPARISON  OF  FIRST  STAGNATION  THEORX  AND  EXPERIMENT 


Experimental  Data  Used 

To  obtain  sufficient  data  for  a  comparison  to  all  the  results  of  the  theory,  it 
was  necessary  to  supplement  the  pressure  and  load  measurements  of  Appendix  H  with 
data  drawn  from  the  schlieren  photographs.  The  information  obtained  from  these  photo¬ 
graphs  was  the  location  and  the  shape  of  the  separation  rocket  jet  boundary  and  the 
floor  stagnation  shock  wave  (see  Figure  8)*  Assuming  the  flow  in  the  expanding  jet 
to  be  correctly  given  by  the  program  of  Reference  (d),  the  flow  conditions  crossing 
d  of  Figure  (8)  at  the  shock  (s)  were  calculated  from  the  photographic  shock  angle. 

Flow  along  the  interstage  floor  was  calculated  from  the  measured  plate  pressure  distri*- 
bution.  Average  flow  conditions  crossing  d,  for  comparison  to  the  theoretical  average 
values,  were  taken  as  the  average  of  the  conditions  at  s  and  f. 

For  the  experimental  configurations  having  simply  the  nozzle  and  the  flat  plate,  the 
jet  boundary  and  the  plate  shock  were  generally  completely  visible  in  the  photographs. 
However,  the  shock  appeared  to  lie  very  close  to  the  plate,  or  at  times  to  disappear 
within  the  plate  on  the  centerline.  This  was  due  to  viewing  the  model  with  a  cone 
of  light  produced  by  the  slit  source  and  returned  by  the  aspheric  mirror  (Figure  3)* 

The  plate  was  thus  viewed  at  an  angle,  and  some  of  the  stagnation  layer  was  obscured 
behind  the  elliptical  plate  shadow.  Although  no  correction  was  required  for  the 
photographic  axial  distances  from  the  nozzle  to  points  on  the  shock,  the  true  position 
of  the  plate  had  to  be  determined  from  the  known  view  angle  and  the  photographic 
position  of  the  rear  edge  of  the  plate.  After  this  position  had  been  determined,  the 
« -Hal  scale  of  the  photograph  was  calculated  from  the  known  and  photographic  separation 
distances. 

For  the  configurations  involving  the  circular  fence,  the  fence  often  obscured  most  of 
the  jet  boundary  and  plate  shock,  except  for  the  view  through  the  region  of  the  axis 
provided  by  a  1.5"  x  2”  plexiglass  window.  The  analysis  of  the  axial  scale  and  hardware 
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positions  vas  as  above,  except  that  both  the  nozzle  and  plate  positions  were  obscured 
and  required  determination  from  the  angles  of  view.  The  areas  of  the  jet  boundary- 
blocked  by  the  fence  were  filled  in  by  noting  the  pressure  to  which  the  known  flow 
in  the  jet  was  raised  in  crossing  the  visible  boundary  shock  wave,  and  continuing 
this  shock  through  the  jet  in  such  a  way  as  to  be  a  smooth  curve  and  maintain  the 
pressure  behind  it  constant.  The  intersection  of  the  plate  shock  with  the  boundary 
was  determined  from  the  derived  boundary  shape  and  the  ring  of  interaction  {repre¬ 
senting  this  intersection  which  was  visible  in  the  window. 

Comparison  of  Theory  and  Experiment 

The  experimental  data  for  three  test  configurations  is  compared  to  the  first  stag¬ 
nation  theory  in  Figures  9-11.  Figures  (9)  and  (10)  refer  to  free  jets,  at  a  nozzle 
pressure  ratio  (Rn/j^)  of  12,000,  impinging  on  flat  plates.  Figure  (11)  is  the  jet 
impingement  within  an  interstage  configuration  at  a  nozzle  pressure  ratio  of  1130. 
All  three  cases  involve  the  sonic  nozzle  of  Figure  (U).  The  solid  lines  are  experi¬ 
mental  values  and  the  dashed  lines  are  from  the  theory  of  Appendix  I. 
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From  these  figures  it  is  seen  that  there  is  excellent  agreement  between  theoretical 
plate  shock  standoff  distance  on  the  jet  centerline  ^ taken  from  Reference  p)  and  the 
standoff  measured  from  the  photographs.  As  a  consequence,  there  is  good  agreement 
between  theoretical  and  experimental  centerline  plate  pressures,  which  are  governed 
by  shock  standoff.  However,  the  theoretical  jet  boundary  (as  predicted  by  the  circular 
arc  method  of  Reference  n)  deviates  greatly  to  the  outside  of  the  photographic  boundary. 
The  radius  of  the  stagnation  layer  exit  (  of  Figure  8)  is  thus  much  greater  in  theory 
than  in  experiment.  The  deviation  of  the  various  theoretical  results  from  experimental 
values  becomes  difficult  to  interpret  in  view  of  the  disparity  in  geometry,  and  the 
fundamental  way  in  which  ^  enters  into  the  analysis  (see  Appendix  I).  It  should  be 
noted  that  theoretical  plate  load  (L),  layer  exit  static  pressure  (  ),  pitot 

pressure  (  1^-  ),  flow  angle  (  oy  ),  and  Mach  number  (  )  agree  remarkably  well 
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with  the  experimental  values  at  the  radius  (  )  taken  by  the  theory,  even 

though  the  actual  layer  exit  is  not  located  there.  However,  the  first  stagnation 
theory  of  Appendix  I  must  be  regarded  as  unsuitable  for  use  on  the  recirculation 
problem  as  long  as  it  fails  to  predict  the  correct  height  and  radius  of  the  stag¬ 
nation  layer* 

Investigation  of  Methods  for  the  Jet  Boundary 

The  method  used  to  obtain  the  jet  boundary  is  not  fundamental  to  the  analysis  of  the 
stagnation  layer  but  an  independent  problem,  since  the  expanding  separation  rocket 
jet  is  supersonic.  The  circular  arc  method  of  Reference  (n)  was  chosen  because  of 
its  speedy  and  because  the  comparisons  to  experiment  in  that  reference  were  favorable* 

'f’he  need  for  speed  in  the  boundary  calculation  for  the  interstage  application  arises 
because  the  cavity  pressure  surrounding  the  jet  is  not  known  a  priori  but  is  an 
iteration  variable.  Each  iteration  necessitates  the  calculation  of  a  new  jet  boundary. 
Tims,  if  the  boundary  calculation  is  time  consuming,  the  length  of  time  necessary  to 
completely  solve  one  Interstage  problem  could  be  prohibitive. 

The  method  of  characteristics  has  been  shown  to  give  boundaries  which  agree  quite  well 
with  experiment  (see  Reference  p).  This  method,  which  produces  very  accurate  results, 
is,  however,  extremely  time  consuming.  An  example  of  the  agreement  between  the  program 
of  Reference  (d)  and  the  present  experimental  data  is  shown  by  the  calculated  points 
in  Figure  (9). 

In  the  interest  of  speed  and  ease  of  calculation,  a  number  of  approximate  methods  have 
been  developed  for  boundary  determinations  including  the  circular  arc  method  of 
Reference  (n)  adopted  in  the  theory.  Three  methods  are  presented  in  Figure  (12),  compared 
to  the  method-of -characteristics  boundary  calculated  for  the  caje  of  Figure  (U).  Good 
agreement  is  seen  for  both  the  Latvala  approximate  method  (also  given  in  Reference  n) 
and  the  Adamson  and  Nicholls  method  (Reference  q)  for  small  values  of  X /rn  (see  inset. 
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Figure  12),  with  the  Latvala  method  being  slightly  better*  Both  methods  are 
more  accurate  in  this  region  than  the  Latvala  circular  arc*  However,  both 
methods  shortly  cross  over  the  characteristics  solution  and  thence  continually 
deviate  in  such  a  manner  as  to  produce  boundaries  considerably  in  excess  of  the 
correct  boundary  and  even  in  excess  of  the  circular  arc  boundary. 

The  circular  arc  approximation  could  be  much  improved  by  decreasing  the  initial 
expansion  angle.  This  decrease  has  no  basis  in  theory,'  but  a  similar  decrease 
was  found  advantageous  by  Latvala  for  similar  pressure  ratios  (Reference  n)*  At 
larger  pressure  ratios,  the  agreement  further  downstream  still  leaves  much  to  be 
desired,  being  generally  larger  than  the  actual  boundary,  then  curbing  too  rapidly 
and  cutting  across  the  boundary* 

It  thus  appears  that  none  of  these  approximate  methods  are  suitable  for  predicting 
the  jet  boundary  a  large  number  of  exit  radii  downstream  of  the  nozzle,  which  is  the 
region  of  interest  in  the  recirculation  calculation.  The  only  reliable  method  is  the 
method  of  characteristics,  which  has  the  large  disadvantage  of  being  oxcessively  time 
consuming.  Either  some  concise  tabulation  of  characteristics  boundaries  for  the  nozzle 
in  question  would  have  to  be  employed,  or  some  new  approximate  method  derived  frota 
characteristics  theory  would  have  to  be  developed,  in  order  to  permit  the  continuation 
of  the  first  stagnation  calculation  and  the  solution  to  recirculation. 
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Investigation  of  the  Assumptions  of  the  Theoi 


To  further  investigate  the  assumptions  of  the  first  stagnation  analysis,  the  theoreti¬ 


cal  jet  boundary  estimates  were  forced  to  agree  with  the  photographic  boundaries  of 


Figures  9-11,  by  adjusting  the  boundary  expansion  angle  at  the  nozzle,  as  suggested 
by  Latvala,  in  Reference  (n)<  The  first-stagnation  analysis  was  then  repeated  for  these 


cases*  The  results  are  shown  in  Figures  (9-11)  by  the  dotted  lines« 
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It  is  Been  that  the  comparison  between  theoretical  and  experimental  stagnation  1  yer 
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geometry  and  flow  properties  is  good  only  for  Figure  (9)*  The  theoretical  layer 
exit  height  (  A Z-d  )  becomes  progressively  greater  than  the  experimental  value 
in  Figures  (10)  and  (11).  T^B  deviation  is  believed  due  to  the  influence  of  the 
jet  boundary.  For  the  case  of  Figure  (9),  the  ambient  pressure  is  low  and  the 
length  of  boundary  involved  is  fairly  small,  so  that  the  boundary  region  (region 
between  jet  boundary  shock  and  the  interface  with  ambient  air)  is  thin,  and  the 
strength  of  flow  in  this  region  is  small.  The  assumption  of  the  theory  that  the 
boundary  is  negligible  appears  good  for  this  case.  For  the  case  of  Figure  (10), 
however,  the  boundary  is  apparent! v  long  enough  to  have  captured  a  substantial  amount 
of  gas  from  the  jet  core.  This  boundary  evidently  prevents  appreciable  flow  from 
entering  the  stagnation  layer  across  the  plate  shock.  Thus  the  theoretical  method, 
which  assumes  all  the  jet  gas  to  be  involved  in  the  layer,  predicts  appreciably  too 
great  a  layer  height.  This  situation  is  carried  to  extremes  in  Figure  (ll),  where  the 
high  ambient  pressure  of  the  interstage  produces  a  very  thick  and  strong  jet  boundary* 

In  fact,  at  the  experimental  layer  exit  point,  sd,  (see  Figure  8)  over  half  the  thrust 
of  the  jet  is  concentrated  in  the  boundary.  Under  these  circumstances,  the  assumption 
that  the  boundary  is  negligible,  and  that  all  the  jet  gas  is  contained  in  the  isentropio 
core  and  crosses  the  plate  shock,  gives  ridiculous  results# 

To  further  investigate  the  theoretical  assumptions,  the  first  stagnation  method  was 
applied  to  the  base  of  Figure  (ll),  with  not  only  the  jet  boundary  adjusted  to  agree 
with  experiment,  but  with  the  axial  momentum  and  mass  flow  considered  to  crws  the 
plate  shock  (C*  and  Cc  of  Appendix  I)  adjusted  to  the  actual  amounts  contained  in 
the  jet  core,  ^he  results  are  shown  in  Figure  (11)  by  the  dot-dash  lines.  The  pre¬ 
dicted  layer  geometry  is  seen  to  be  very  good,  as  is  layer  exit  flow  angle.  However, 
the  plate  load  is  considerably  lower  in  theory  than  in  experiment,  the  theoretical  layer 
exit  Mach  number  (  Hd  )  Is  much  too  high  and  the  exit  static  pressure  (  )  is  much 

too  low*  Since  these  results  were  not  obtained  in  the  adjusted-boundary  calculation  of 
Figure  (9),  they  point  to  an  influence  of  the  relatively  high-density  -  oundary  in 
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restraining  the  exit  of  flow  from  the  layer,  causing  a  build-up  of  pressure  within 
the  layer  and  on  the  plate*  It  appears  that  the  calculation  of  first  stagnation 
of  the  Jet  on  the  interstage  floor  cannot  be  made  by  ignoring  the  jet  boundary  and 
assuming  the  jet  flow  to  be  isentrcpic.  In  the  high  pressure  interstage  environment, 
the  requirement  that  the  high-density  boundary  stream  be  turned  by  the  flow  accelerating 
across  the  plate  out  of  the  stagnation  layer,  evidently  imposes  a  different  character 
to  the  layer  from  that  to  be  expected  in  near-vacuum  impingement*  ^he  presence  of  the 
jet  boundary  must  be  incorporated  into  the  first  stagnation  theory  before  good  results 
can  be  expected* 

It  should  be  noted  from  Figure  (11 )  that  the  pressure  gradient  across  the  interstage 

floor  is  sufficiently  great  to  induce  supersonic  flow  between  the  first  stagnation 

layer  exit  and  the  sides  of  the  interstage.  Thus,  the  first  stagnation  process  is 

independent  of  the  side  stagnation  process.  Likewise,  supersonic  flow  exists  along 

the  sides  between  the  bottom  and  top  of  the  interstage.  It  would  seem  that  the 

assumption  of  a  step-wise  sequence  of  independent  stagnations  is  correct,  and  that 

a  theoretical  method  for  recirculation  which  is  based  on  this  assumption  would  be 
ic 

feasible. 
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The  question  of  how  the  high  density  jet  boundary  of  Figure  (11 )  (total  pressure  - 
260  psi)  interacts  and  is  turned  along  the  plate  by  the  relatively  low  density  stag¬ 
nation  layer  flow  (total  pressure  «  U  psi)  is  fundamental  to  an  understanding  of  the 
conditions  which  would  have  to  be  incorporated  into  the  first  stagnation  theory  to 
make  it  correct.  This  question  has  studied  by  a  conservation  analysis  of  the 

small  volume  shown  in  the  inert.,  Figurt  (11),  making  use  of  the  measured  plate  pressure 
distribution  and  assumed  linear  pressure  and  flow  angle  distributions  across  the  flows. 
In  this  analysis,  the  boundary  and  layer  flows  have  been  assumed  not  to  mix,  but  to 
rotain  their  separate  identities.  Since  the  layer  flow  is  subsonic,  it  has  been  assumed 
isentropic.  The  results,  shown  in  Figure  (11),  are  that  the  boundary  flow  turns,  and 
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grows  greatly  in  width,  while  losing  total  pressure  rapidly.  The  layer  flow  con¬ 
tinues  to  accelerate  towards  M  “  1,  and  thus  becomes  smaller  in  thickness.  The  indi¬ 
cation  is  that  as  flow  continues  to  expand  to  greater  radii,  the  boundary  stream  will 
continue  to  broaden,  turn,  and  lose  total  pressure,  occupying  an  ever  increasing  pro¬ 
portion  of  the  ever  decreasing  total  layer  width.  The  method  of  turning  the  boundary 
must  be  through  shocks  induced  in  that  stream  by  the  subsonic  layer  flow  underneath, 
and  reflected  expansions  induced  when  these  shocks  strike  the  interface  with  ambient 
air*  This  process  continues  until  the  layer  flow  reaches  M  *  1,  at  which  point  shocks 
can  propagate  al3o  through  this  stream  and  reflect  from  the  plate.;  The  hump  in  the 
plate  pressure  distribution  of  Figure  (11 )  is  very  suggestive  of  such  a  process.  The 
result  is  expected  to  be  a  rather  thin,  two -layered  high  density  flow  along  the  plate, 
with  the  boundary  flow  occupying  most  of  the  width. 

One  further  observation  on  the  theoretical  assumptions  arises  from  an  attempt  to  apply 
the  first  stagnation  theory  to  near-vacuum  impingement  on  the  cone-shaped  experimental 
plate  configuration.  The  method  or  Reference  (o)  predicted  a  centerline  shock  standoff 
distance  much  smaller  than  the  experimental  value.  Using  this  distance,  the  calcu¬ 
lations  beccsne  mathematically  unstable  and  no  solution  was  found.  The  instabilities 
were  of  a  type  previously  encountered  in  the  development  of  the  theory  (see  Appendix  I), 
and  are  believed  to  be  associated  with  insufficient  stagnation  layer  cross-sectional 
area  during  iteration.  With  the  small  centerline  standoff  distance,  the  theory  could 
not  pass  a  spherical  shock  out  the  cone  and  achieve  sufficient  layer  area  at  all  points® 
This  is  shown  in  Figure  (lj).  Thus,  the  method  of  Reference  (o)  is  not  adequate  to  give 
centerline  standoff  distance  when  the  interstage  floor  has  a  slope  at  the  center.  This 
is  not  surprising  since  the  assumpti  _>r.  that  the  plate  surface  near  the  jet  centerline  is 
normal  to  the  centerline  is  fundamental  to  the  analysis  in  Reference  (o).  That  method 
would  have  to  be  generalized  for  sloped  surfaces  in  order  to  handle  such  interstage  floor 
shapes.  Thi3  conclusion  may  also  have  important  bearing  on  jet  stagnation  calculations 
when  the  jet  is  tilted  with  respect  to  the  interstage  floor  axis® 


IMPLICATIONS  OF  THE  THEORETICAL-EXPERIMENTAL  COMPARISONS 


The  important  implication  to  design  studies  of  the  experimental  and  theoretical  re¬ 
sults  is  that  in  rocket  staging  problems  the  jet  gas  from  the  separation  rocket  may 
be  confined  to  a  thin,  high  density  stream  as  it  flows  along  the  interstage  floor* 
Thus,  by  supplying  vents  in  the  sides  of  the  interstage  at  the  level  of  the  floor, 
much  of  the  jet  gas  can  be  exhausted,  and  the  interstage  recirculation  (and  hence 
staging  forces)  can  be  greatly  reduced*  However,  if  vents  are  above  the  floor,  even 
by  a  small  amount,  they  may  be  much  ies3  effective,  because  the  thin  floor  stream  may 
be  expected  to  require  only  a  small  height  in  which  to  complete  its  turn  at  the  sides 
and  begin  re-accie rating  toward  the  top  of  the  interstage*  For  a  design  having 
simultaneous  rocket  ejection  of  two  or  more  adjacent  stages,  strong  flows  exiting 
from  any  bottom  vents  would  interact  to  produce  a  potentially  strong  flow  up  be¬ 
tween  the  stages,  which  might  produce  high  staging  perturbations.  This,  of  course, 
is  the  dase  for  the  Polaris  A3  (Reference  i)«  Under  these  circumstances  It  might  be 
better  to  avoid  bottom  venting  entirely* 
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CONCLUSIONS 


1«  A  ten  week  test  program  including  2 59  successful  runs  was  conducted  at  the 
LMSC  Tenney  Altitude  Facility.  High  pressure  nitrogen  was  used  to  simulate 
flow  from  a  rocket  motor  exhausting  into  the  test  cell  at  pressure  altitudes 
which  varied  from  105*000  feet  to  117*000  feet.  Testing  was  completed  on 
all  the  configurations  planned  with  adequate  pressure  surveys,  loads,  and 
photographic  data  obtained  to  permit  theoretical  analysis* 

2*  The  recirculation  flow  within  an  Interstage  during  rocket  powered  stage 
separation  appears  to  consist  of  a  chain  of  supersonic  accelerations  and 
stagnations,  each  of  which  may  be  solved  independent  of  the  others.  A 
method  appears  feasible  for  solving  staging  problems  involving  iteration 
on  interstage  cavity  pressure,  with  calculation  of  the  chain  of  stagnations 
around  the  interstage  for  each  iteration. 

3.  The  first  stagnation  of  the  separation  rocket  jet  on  the  interstage  floor 

appears  to  follow  the  laws  of  conservation  of  mass,  axial  momentum  and  radial 
momentum,  as  applied  to  the  entire  stagnation  layer.  A  method  for  calculating 
this  process  with  these  equations  appears  feasible  providing: 

a.  The  correct  position  of  the  jet  boundary  is  known; 

b.  The  mass  and  momentum  in  the  jet  is  properly  apportioned 
between  the  isentropic  core  and  the  boundary  strip; 

c.  The  influence  of  the  jet  boundary  in  restricting  the  flow 
across  the  exit  of  the  stagnation  layer  is  properly  considered. 

he  In  recirculation  cases,  the  gas  flow  across  the  interstage  floor  appears  confined 
to  a  thin  layer  due  to  the  action  of  the  jet  boundary  on  the  first  stagnation  flow 
field.  This  would  be  expected  to  render  interstage  vents  at  the  level  of  the 
floor  very  effective  in  reducing  recirculation  forces.  Likewise,  for  a  cluster 
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of  separating  stages,  sealing  each  interstage  at  the  level  of  the  floor 
would  appear  to  be  very  effective  in  preventing  recirculation  in  the  center 
of  the  oluster* 


w 
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APPENDIX  I 


CONTROL  VOLUME  METHOD  FOR  AN  AXLSTMMETRIC  FIRST  STAGNATION 
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Basic  Assumptions 

As  shown  in  Figure  (8),  an  axlsymmetric  free  jet  is  assumed  to  expand  from  the 
separating  body  nozzle  to  the  interstage  pressure  environment,  and  stagnate  against 
an  axlsymmetric  interstage  floor  having  the  same  axis  as  the  nozzle  and  jet.  The 
gas  is  assumed  perfect,  with  ratio  of  specific  heats  (  7f  )  constant  at  the  nozzle 
chamber  value.  The  entire  flow  field  is  assumed  steady.  The  jet  expansion  and 
stagnation  processes  are  assumed  adiabatic,  with  the  flow  in  the  jet  isentropic 
except  for  a  small  region  near  the  boundary-.  The  flow  properties  in  the  free  jet 
can  thus  be  obtained  with  the  method  of  characteristics  (Reference  d).  The  flow 
in  the  layer  (V)  is  also  assumed  isentropic,  so  that  all  entropy  changes  in  the 
field  occur  at  the  shock,  8. 

Basic  Equations 

The  basic  equations  used  are  conservation  of  the  mass,  axial  (x)  momentum  and 
radial  (r)  momentum  of  the  gas  within  control  volumes.  These  equations  hold  if  the 
flow  is  steady,  adiabatic,  and  perfect.  Under  these  conditions  the  conservation  of 
energy  is  not  independent  of  conservation  of  mass  and  momentum.  The  mass  and  momentum 
equations  are  written  for  Volume  V  of  Figure  8,  with  mass  and  axial  momentum  fluxes 
into  V  determined  from  the  fluxes  out  of  V-^. 

Since  the  flow  is  assumed  steady,  conservation  of  mass  in  V  and  V-^  takes  the  form* 
(mass  flow  into  volume)  ■  (mass  flow  out  of  volume) 

For  V^  (see  list  of  symbols): 

_  QPt„ 

For  V,  rbfn  =  cut  from  V^.  Thus: 

where  p^,  ,  and  Md  refer  to  average  values  across  the  stagnation  layer  at  d 


z 


bum 


(as  defined  in  the  list  of  symbols),  and  $  has  been  divided  out# 


Conservation  of  axial  momentum  in  the  volumes  has  the  form* 

(flux  of  x -momentum  into  volume)  +  (integral  of  all  pressure 
forces  on  the  volume  in  the  positive  x-direction)  - 
(flux  of  x-tu  omentum  out  of  volume)  “0  , 

since  the  flow  is  steady  and  no  external  forces  are  applied  except  along  the 
bounding  surfaces  of  the  volumes*  for  Volume  Vt  t 

(flux  of  x-momentum  out  of  volume  +  pressure  force  on  shook  surface,  a, 
in  the  negative  x-direction)  ■ 


9  ( /  +  T  Mf)  Pm 

z(l  +  *■=} 


■h 


s. 

z 


Pa 


For  volume  V,  the  flux  of  x-momentun  in  over  s  +  the  pressure  force  in  the 
positive  x-direction  on  a,  is  equal  to  the  left-hand  side  of  the  above  equation* 
Hence: 


2* 


( I+rMn z)  Pr„  nfi 
(/v  Vg  M„z  )  A 


+■  netf-rf-) 


—  Z  Y  Pj  Ax j  sin  cKj  cos  Mj z  =  O 


where  9/2  has  been  divided  out*  Here,  J  Pp  rdr  is  the  pressure  force  exerted  on 
V  in  the  negative  x-direction  by  the  interstage  floor* 


Under  the  basic  assumptions,  conservation  of  radial  momentum  in  volume  V  has  the 

f0«3! 

(flux  of  r-m omentum  into  V  across  surface  s) 

+  (positive  r  component  of  the  pressure  force  acting  on  s) 

+  (positive  if  component  of  the  pressure  force  on  the  sides,  m,  of 
7)  +  (positive  r  component  of  the  pressure  force  exerted  by  the  inter¬ 
stage  floor,  f )  -  (flux  of  r-momentum  out  of  V  across  d)  -  (r  component  of 
pressure  force  acting  on  the  exit,  .d,  of  V)  *  0 
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Denoting  the  sun  of  the  first  two  terms  by  Ijr-  ,  we  have: 

3*  Tjr  +  Z  S'fff  I  Pf  (x>r)  dxJr  ^ 

*4 

-  Pcj  A*c/  e  ('  ^  ? Wj2sinz  <*4) 


e  p£  <r)  dx 

*fc 

=  O 


Equations  1,  2,  and  3  are  the  basic  equations  used  in  the  solution  to  the  first 
stagnation  problem.  The  desired  answers  are  the  flow  properties  and  geometry 
(  J7j  ;  ^ J  xj  )  ihe  exit,  d,  of  the  stagnation  layer, 

and  the  total  load  on  the  interstage  floor  under  the  layer  (/ ZfT P^.{f~)rdr)9 


Assumptions 

In  addition  to  the  desired  flow  properties  (  ,  M  j )  cX^  ),  the  geometry  of 

the  stagnation  layer  exit  (  ty  j  Axj)i  and  the  interstage  floor  pressure  distri¬ 
bution  (  Pj:  (r)  ),  equations  1-3  contain  the  radial  thrust  of  gas  crossing 

s  (Tjr),  the  static  pressure  distribution  in  the  i-yer  (  Pj  ( X,  r)  ),  and  the  shock 
shape  (xg(r))»  There  are  thus  9  unknown  quantities  and  functions  to  be  determined 
with  3  equations.  This  requires  assumptions  permitting  6  of  the  unknown  quantities 
to  be  calculated  or  related  to  other  unknowns* 


If  the  jet  boundary  shape  is  known,  r^  is  related  to  A  Xj  by  simple  geometry  in¬ 
volving  the  known  interstage  floor  shape.  In  this  method  the  jot  boundary  is  assumed 
to  be  given  by  a  Latvala  circular  arc  (Reference  n),  as  follows  (see  Figure  8  and 
list  of  symbols): 


~  (My  )  ~  P  (My,)  -h  cX„ 


where: 


and: 


V  (M) 


=  VW  +°"~w %  w-~) 


'  -  'I 


r 

h  iMu.  -  . 


+  13.78 

372 — - IA 


{r 4-1 )  (  ST  v-  ) 


It 


In  order  to  relate  the  shock  shape,  Xs  (r)  }  to  AXj  and  Qj  >  i*  is  necessary 

to  know  the  point,  XSC  ,  where  the  shock  crosses  the  axis  of  symmetry,  and  also 
the  shape  function  of  the  shock.  In  this  method,  /  is  assumed  to  be  given  by 
the  loshihara  analysis  (Reference  o): 


r  r  _  *se  r  r-z 
6*  “  T~  L  r+l 


'SC 


~  *fc  ~  ^ 


•h 


) 
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where  Mjc(  X sc  )  is  "the  Mach  number  in  the  free  jet  at  the  point  where  the  shock 
crosses  the  axis,  as  given  by  the  method-of  —characteristics  program  of  Reference  d» 
Since  Mj  varies  with  X sc  >  aquations  6  must  be  solved  iteratively. 


id 


The  shape  of  the  shock  is  assumed  in  this  analysis  to  be  a  spherical  cap,  with 
center  on  the  axis  of  symmetry.  This  assumption,  together  with  the  locations  of 

two  points  on  the  shock  (  X sc  9  and  *sd  9  )*  is  sufficient  to  determine  Xs  (  r)t 

-  (*sc  -  *  SJ  )z  -/•  Q Z 


-/ 


P  _ 

/5  Z  (*sc  -  Xsct) 

(f^  (r)  —  S'  in 

Xs  (r)  -  f>s  L^oS  />s(r)  -  /] 


fs 


sc 


Where  <fis  is  the  angle  between  the  normal  to  the  shock  at  the  point  (  X5  }  r  ),  and 


the  axis  of  symmetry* 


With  the  shock  shape  known,  Tjr  may  be  determined,  since  the  flow  properties  in 


the  free  jet  are  given  by  method-of' -characteristics : 

fis  cl 

Tjr  =  f*  e(  Pj  i  ?Mfcos  sin  £  sin  dj  +  ■  sin ]J  <f>s 


8, 


here,  Pj,  °C j  ,  and  Mj  are  the  flow  properties  in  the  free  jet  at  the  point 


X  =  *SC  i  f>  LCOS  tf^  -/  J  > 


r  - 


fs 


s  m 
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Since  changes  in  entropy  (and  hence  total  pressure)  are  assumed  to  occur  only 
at  the  shock,  s,  a  knowledge  of  the  shock  shape  and  all  upstream  jet  flow 
properties  also  pennits  calculating  the  total  pressure  of  flow  after  crossing 
the  shock  at  each  point  on  the  shock: 

?  / 

p  _  p  f  (#+/)  Mj^cos2-  lr-t  f _ r-h  / _  Jr-/ 

Ts  ^LCr-o  Mjco^  (c/j  -  ps  )  4-  z  ]  L2VMj*-co^(cij-<i>s)  -r  +  /  J 

Here,  pj,  Mj,  and  c/j  are  jet  flow  properties  given  as  above.  These  data  permit 
the  calculation  of  the  mean  total  pressure  of  the  layer  flow,  ,  by  scma  sort 

of  averaging  technique.  In  this  method  an  area -weighted  average  is  taken* 

9. 

o 

Using  Prd  ,  Md  and  Pd  may  be  related  by  the  equation: 


J  R~s  cf>s  cos  (cfj  -  dcf>5 
J  Sei  Sin  (dj  -  <4  )  d  <ps 


Finally,  the  mean  flow  angle  at  the  layer  exit,  <*j  ,  may  be  related  to  the  shock 
shape  and  ixy  if  the  shape  of  the  flow  angle  distribution  across  the  layer  is 
assumed.  This  is  because  the  shock  shape,  together  with  the  known  jet  flow  field 
upstream  specifies  the  angle  of  flow  just  downstream  of  the  shock  at  all  points  along 
the  shock,  including  at  the  layer  exit: 


Z  +«n  (djj ~  <j>sd)  LMj/cos-  (otj’j  -  <f)iet)  -  /  J 
Z  4-  Mjf  1*  +  l  -  Z  cos*  (*jj  -  4sd  )  J 


(Here  c/jj  ,  Mjd  are  jet  flow  properties  just  upstream  of  the  shock  at  the  intersection 
of  the  shook  with  the  jet  boundary).  Algo,  U>e  f.1^;  a^igle  at  the  layer  exit  on  the 


interstage  floor,  dr, ,  must  fcn  the  same  as  the  angle  of  the  floor  itself*  Thus,  the 
flow  angle  £ '*  the  two  e*?ds  of  the  segment  A  Xj  is  specified  by  the  shock  shape*  In 
this  analysis  the  flow  angle  distribution  across  the  layer  at  d  is  assumed  linear# 
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Thus* 

11.  Ctj  =  ±  (C<SJ  +  ctfj  ) 


In  this  method  the  interstage  floor  pressure  distribution  is  replaced  by  an  average 
floor  pressure,  Pf  ,  which  gives  the  load  on  the  interstage  floor: 

rj 

12.  j  ?f(r)rdr  =  ^  rj2- 

Xf  /  /Xfci 

(  p,cr)  r(x)  dx  -■  If  J  (x) 

JHc  *fc 

Furthermore,  the  layer  static  pressure  distribution  is  replaced  by  an  average  static 
pressure,  which  is  assumed  equal  to  Pf : 

(rJ  fXf  r  rcl 

13*  j  J  T>U,r)JxJr  =  fyj(Xf-Xs )  Jr 


Equations  U  -  13  reduce  the  number  of  basic  unknowns  in  equations  1  -  3  to  .  A  XJ , 
Md,  and  Pf*  Therefore,  there  are  sufficient  basic  equations  to  determine  these 
remaining  unknowns.  However,  the  equations  are  very  difficult,  due  to  the  term 
(  sin  ^  ).  Since  G  is  arbitrary,  it  has  been  convenient  to  assume  it  to  be  a 
small  value,  so  that  sin  S  ~  *  When  this  is  done,  Q  divides  out  of  equation  3» 

and  this  variable  is  eliminated  frcm  the  entire  problem. 


Method  of  Solution 


In  this  analysis  the  boundary  shape  is  determined  from  equations  Ij  and  5,  and  the 
shock  centerline  standoff  is  determined  from  6,  using  the  given  nozzle  exit  conditions, 
ambient  pressure,  and  interstage  floor  shape*  The  procedure  is,  next,  to  enter  a 
Axj  iteration.  A  XJ  is  assumed,,  and  from  this  r^,  the  shock  shape,  Tjr,  Pj^, 

rJ  .  rKfd 


(Xf  -  X$  )  cfr  and J  Pp  dx  »  are  determined 

V  A 

equations  1  -  3>  together  with  10,  are  solved  for  pf,  Md,  and  a  new  l\Xj  •  Ihe 
iteration  is  continued  until  the  assumed  and  calculated  A/J  agree.  _ 


Then. 
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It  has  been  found  that  this  procedure  is  mathematically  very  unstable  when 

equations  1  -3  are  used  without  modification.  For  setae  assumptions  of  A  XJ 

there  is  no  solution  for  Md;  for  other  assumptions  a  value  for  Md  is  obtained, 

but  the  calculated  Ax^may  be  negative  or  imaginary.  This  behavior  has  been 

r  rj 

traced  to  the  term  t  ^  y  (  Xy  -  X  s  )  ci  r~ .  For  a  range 

of  assumptions  of  Aty  which  are  too  small,  this  term  is  too  small  to  satisfy 
equation  3  with  flow  being  out  of  the  stagnation  layer  at  d,  as  implied  by  equations 
1  and  11.  It  has,  therefore,  been  found  convenient  to  increase  the  size  of  this 
term  in  case  the  assumed  Axj  is  too  small,  so  that  the  equations  become  mathe¬ 
matically  stable.  Consequently,  the  term: 

Arn  c  ^ ' a  Siam  eo/  ^^^ho  be  cot  !c.u  la  lec/ 

(rJ 

has  been  added  to  Pf J  (Xr  ~X  $)  a  r  in  equation  3*  This  terra  compensates  for 

t )  ' 

too  small  a  value  of  (assumed),  yet  disappears  when  AXj(  as  sinned)  approaches 

the  calculated  value.  The  form  of  this  term  arises  from  a  consideration  of  the 
effect  of  an  error  in  hxj  on  the  layer  area,  J  fxr  -  Xs)  dr  . 

Eliminating  Pd,  Pf,  and  A  Xj  (  calculated)  from  equations  1,  2,  10,  and  the  modified 
equation  3j  we  have: 

^  Pd  =  V * 


asuimcc/ 


&XJ 


■ho  be  cd  !c.u  la  le.d 


Cc.  K  /  o  *  _  2vCi  COSQ/J  Me!  v 


Within  each  &XJ  iteration,  equation  XU  is  used  to  solve  for  Md.  Pd  as  calculated 
from  both  expressions  of  lii  is  forced  to  be  the  same  by  trial  and  error  adjustment 
of  Md.  When  a  value  of  Md  for  the  iteration  lias  been  found,  a  calculated  hxj  is 
obtained  from  equation  It 


1. 


Cl 


%  ry  sinclj  Mj*-  ' 

When  the  &XJ  iteration  has  converged,  Pf  is  calculated  by: 

_  -//  +  lOt  +■  ?a  m*-  )  -  2fCt,  coscQ 

V*  //  y-  MJ-' 

and  the  total  load  on  the  interstage 


?  = 


floor  is  given  by: 
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Explanation 

The  data  presented  herein  is  divided  into  two  broad  classifications,  sonic  nozzle 
configurations  and  supersonic  nozzle  configurations.  The  pressure  data  is  normalized 
to  a  nozzle  total  pressure  of  1000  psia.  To  make  use  of  the  pressure  data  for  a 
different  total  pressure,  multiply  the  pressure  data  by  the  new  total  pressure  times 
10~3.  The  force  ratios  can  be  used  directly. 

All  the  data  presented  was  taken  at  the  same  value  of  nozzle  total  pressure  to  ambient 
pressure  ratio,  1.2  x  10^*  For  "closed"  configurations,  that  is,  configurations  having 
a  nozzle  plate  and/or  a  circular  fence,  the  local  ambient  pressure  is  not  the  same  as 
the  ambient,  or  test  cell,  pressure.  In  particular,  for  configurations  having  a 
nozzle  plate,  the  local  ambient  pressure  is,  in  general,  vastly  different  from  the 
test  cell  pressure.  This  local  ambient  pressure,  which  determines  the  shape  and 
extent  of  the  jet  boundary,  is  essentially  unknown.  Hence,  the  pressure  ratio 
mentioned  above  is  based  on  the  test  cell  pressure. 

The  circular  fence  used  throughout  the  test  had  a  height  of  21»7  exit  radii  for  the 
< 

sonic  nozzle,  and  9*07  exit  radii  for  the  supersonic  nozzle.  The  inside  diameter 
of  the  fence,  the  base  plate  diameter  and  the  nozzle  plate  diameter  were  all  86.8 
exit  radii  for  the  sonic  nozzle,  and  36.28  exit  radii  for  the  supersonic  nozzle. 

The  locations  of  the  pressure  probes  used  in  the  test  are  indicated  on  the  plots  of 
rake  pressures.  The  symbols  H,  Rg,  and  ^  ,  used  to  locate  the  probes,  are  defined 
as  in  Figure  (6),  as  are  all  other  symbols  encountered  in  the  plots# 

Summary  of  the  Data 


fr„  *  1000  psia 

Supersonic  Nozzle: 

V* 

-  12000  Y  -  1.1* 

a/a*  * 

<<«  -  16° 

Sonic  Nozzles 

A/A*  *  1#0  -  Oo 
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Figure  1 


Figure  2 


Figure  3 


Figure  U 


Figure  !? 


Sonic  Nozzle  Flat  Plate 

a#  Plate  Pressure  Distribution 

b.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

Sonic  Nozzle  Flat  Plato  with  Nozzle  Plate 

a*  Plate  Pressure  Distribution 

b.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

c.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation 
Distance  and  Nozzle  Plate  Pressure  Distribution 

Sonic  Nozzle  Flat  Plate  with  Circular  Fence  (h /rn  m0) 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

Sonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and 

Circular  Fence  (  hjrn  -  0) 

a*  Plate  Pressure  Distribution 

b*  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

d.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Distance 
and  Nozzle  Plate  Pressure  Distribution 

Sonic  Nozzle  Flat  Plate  with  Circular  Fence  ( h/r„  -10190) 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs  Separation 
Distance 


Figure  6 


Figure  7 


Figure  8 


Figure  9 


Sonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and 

Circular  Fence  (  h  j r„  »  10.98) 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation  Distance 

d.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation 
Distance  and  Nozzle  Plate  Pressure  Distribution 

Sonic  Nozzle  •  Flat  Plate  with  Circular  Fence  (h/r„  «l6.30 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation  Distance 

Sonic  Nozzle  Flat  Plate  with  Circular  Fence  (  hjr„  -21.70 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  rs 

Separation  Distance 

Sonic  Nozzle  Simulated  Re-entry  Body  Flare 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation  Distance 

d.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Distance  and 
Nozzle  Plate  Pressure  Distribution 
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Figure  10 


Figure  11 


Figure  12 


Figure  13 


Figure  11* 


Sonic  Nozzle  Flat  Plate  with  Nozzle  Tilted  U»?°  aud  9*8° 

a.  Plate  Pressure  Distribution  -  Plane  1 

b.  Plate  Pressure  Distribution  -  Plane  2 

c.  Plate  Pressure  Distribution  -  Plane  3 

d.  Plate  Force/Thrust  Ratio  and  Rake  Pressures 

vs  Separation  Distance  -  U*7U 

e«  Plate  Force  Thrust  Ratio  and  Rake  Pressures 
vs.  Separation  Distance  -  9*8° 

Sonic  Nozzle  Cone,  Altitude  ■  21*7  radii 

a.  Plate  Pressure  Distribution 

b.  .  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation  Distance 

Supersonic  Nozzle  Flat  Plate 

a*  Plate  Pressure  Distribution 

b*  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate 

a*  Plate  Pressure  Distribution 

b«  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

c*  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Distance 
and  Nozzle  Plate  Pressure  Distribution 

Supersonic  Nozzle  Plat  Plate  with  Circular  Fence  {^/rn  *  0) 

a.  Plate  Pressure  Distribution 
b*  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 


Figure  V~> 


Figure  16 


Figure  17 


Figure  18 


* 


Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and 

Circular  Fence  (  h/rn  *0) 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

o*  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation 

d.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Distance 
and  Nozzle  Plate  Pressure  Distribution 


Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and 

Circular  Pence  (  h  Jrn  ■  1.13U) 

a.  Plate  Pressure  Distribution 

b*  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressure 

vs  Separation  Distance 

d.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation 
Distance  and  Nozzle  Plate  Pressure  Distribution 
« 

Supersonic  Nozzle  Flat  Plate  with  Circular  Fence  (^/r„  “2.268' 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation  Distance 


Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and  Circular 

Fence  (  h/r„  -  2.268) 

a.  Plate  Pressure  Distribution 

b.  Fence  Pressure  Distribution 

c.  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 

Separation  Distance 

d.  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Distance  and 
Nozzle  Plate  Pressure  Distribution 
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Figure  19 


Figure  20 


Figure  21 


Figure  22 


Supersonic  Nozzle  Flat  Plate  with  Circular  Fence  (  h/i~n 
a*  Plate  Pressure  Distribution 
b*  Fence  Pressure  Distribution 

l 

c*  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
.‘.Separation  Distance 

Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and  Circular 

Fence  (  h/r*  -  U*535)  j 

a*  Plate  Pressure  Distribution 

b*  Fence  Pressure  Distribution 

o*  Plate  Force/lhrust  Ratio  and  Rake  Pressures  vs  ' 

Separation  Distance 

d*  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Distance 
and  Nozzle  Plate  Pressure  Distribution 

Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and  Nozzle 

Tilted  10° 

a.  Plate  Pressure  Distribution  -  Plane  1 

b.  Plate  Pressure  Distribution  -  Plane  2 

c.  Plate  Pressure  Distribution  -  Plane  3 

d*  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

e*  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation 
Distance  and  Nozzle  Plate  Pressure  Distribution 

Supersonic  Nozzle  Flat  Plate  with  Nozzle  Plate  and  Circular 

Fence  -  0)  Nozzle  Tilted  10° 

a*  Plate  Pressure  Distribution  -  Plane  1 

b.  Plate  Pressure  Distribution  -  Plane  2 

c.  Plate  Pressure  Distribution  -  Plane  3 

d«  Fence  Pressure  Distribution 

@«  Plate  Force/Thrust  Ratio  and  Rake  Pressures  vs 
Separation  Distance 

f*  Nozzle  Plate  Force/Thrust  Ratio  vs  Separation  Dis¬ 
tance  and  Nozzle  Plate  Pressure  Distribution 


